1, 2008 5:11 WSPC/INSTRUCTION FILE Beijingpenta 



International Journal of Modern Physics A 
© World Scientific Publishing Company 



A Review of Pentaquark Calculations on the Lattice 



Kch-Fei Liu 

Dept. of Physics and Astronomy, University of Kentucky, Lexington, KY, 40506, USA 

Nilmani Mathur 

Jefferson Lab., 12000 Jefferson Ave., Newport News, VA 23606, USA 



We review lattice calculations of pentaquarks and discuss issues pertaining to inter- 
polation fields, distinguishing the signal of pentaquarks from those of the KN scattering 
states, chiral symmetry, and ghost state contaminations. 

Keywords: Pentaquarks; lattice calculation 



1. Introduction 

The recent interest in pentaquark baryons has been inspired by the experimental 
discovery of O + (1540) whose quark composition is uudds. We will not address the 
experimental situation. It is summarized by T. Nakano in his plenary talk during 
this conference 0. 

In the past few years, there has been quite an effort to calculate these pentaquark 
states to see if they can be observed in lattice QCD calculations and if their masses 
and other properties agree with those observed experimentally. In particular, it 
would be interesting to see if there is a way to understand why the observed width 
of + (154O) is much narrower than those of the ordinary baryons which are of the 
order of several hundred MeV. 

Somewhat similar to the experimental situation, among the dozen or so lattice 
calculations, 5 calculations claim to have positive signals for the pen- 

taquark; while 7 calculations l 4 l v l v l ffl i v l lz l i '.? l reported null results. Prior to pen- 
taquarks, lattice calculations dwell mostly on stable particles in strong decays, such 
as the nucleon, the pion, and kaon. Since the quark mass in the present lattice calcu- 
lations is still not light enough, the A, p, and <f> are also below the decay thresholds 
and can be calculated as the ground states. Come pentaquark, the interpolation field 
involves two u quarks, two d quarks, and one strange antiquark which inevitably will 
couple to the KN scattering states in addition to the potential pentaquark state. 
This presents a challenge not confronted before and, presumably, posted some con- 
fusion in the first round of lattice calculations. We shall address some of the issues 
due to this complication, namely the questions of interpolation fields, the issue of 
KN scattering states, the importance of chiral symmetry, and the contamination 
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of ghost states - a quenched artifact and will attempt to draw some lessons from 
these calculations. 

2. Interpolation Fields 

There has been quite a lot of discussion in the literature on the choice of inter- 
polation fields. The most naive one is the product of nucleoli and kaon interpolation 
fields, i.e. 

x /=0,l = e abc (^a^^) (-e^e ) ^ {u „ d}] > (1) 

where the =p combination is for 7 = and 1=1 and has been denoted as such.. 
Another one is similar, except with the color indices e and c positioned differently 

x 7=0,l = £ abc (^Ta^^) [y e (-e^e ) ^ {u „ d}] ^ (2) 

The third one which has been used in the lattice pentaquark calculation is 

x r = e g C e eg fh e abc ^T.^) ( U Tf CTd h ) TC^S^, (3) 

with r = {S, A} = {1,7^75} ^. It is motivated by the diquark-diquark-antiquark 
picture of Jaffe and Wilczek 

The masses and spectral weights due to an interpolation field operator O can be 
extracted from the zero-momentum correlation function with a point source, which 
is a sum of exponentials for the spectrum of states with the quantum numbers of 
the operator O 

{J2O(S,t)O(S ,0))=J2Wie- mit . (4) 

x i 

It is often asserted in the literature that xij being the KN interpolation field, will 
couple to the KN scattering states stronger than do \i an d X3 ■ O n the other hand, 
X3 is expected to couple weakly to KN states and more strongly to the pentaquark 
state. However, it has been learned in the lattice community over the years not 
to be overly sanguine about predicting the structure of the hadron based on the 
interpolation field. It is already shown in Ref. [4] that Xi,X2, and X3 are linearly 
related through the multiplication of 75 and a Fiertz transform between the u and 
s fields, i.e. 

75X( X (=°- X r )^(xf + Xi)- (5) 

Thus they are expected to couple to the same states (1/2 - and l/2 + ) with compa- 
rable spectral weights. 

As a further check, we show in Fig. 1 the effective mass plot of the ground 
states as obtained with the three operators in the l/2~ channel for the u/d quark 
mass corresponding to a pion mass at 633 MeV and the strange mass corresponding 
to the physical (j> mass on a 16 3 x 28 lattice with overlap fermions and a lattice 
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Fig. 1. The effective mass of the 1/2 ground states from the Xi>X2iX3 an d the KN open jaw 
diagram as a function of time separation. 
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Fig. 2. The spectral weights (in arbitrary units) of the 1/2 (left panel) and l/2 + (right panel) 
ground states from the XiiX2i and X'i interpolation fields. 



spacing at 0.2 fmEl We see that at short time separation, i.e. t < 7(1.4fm), the 
three interpolation fields do not give the same mass. Only when t ~ 12 (2.4fm) do 
they come down to the KN threshold. Whereas the 'open-jaw diagram', where the 
kaon and nucleon are separately projected to the zero momentum states, has better 
overlap with the threshold NK scattering state and come down to the threshold 
much earlier. This shows that the three interpolation fields project to the same 
ground state at large time separation. We also plot their ground state spectral 
weights in the 1/2 - and the 1/2+ channels in Fig. 2. We see that the spectral 
weight of X3 is somewhat smaller than those of xi an d X m the l/2~ channel 
{KN in S-wave), but not orders of magnitude smaller; whereas, all three give the 
same weights in the 1/2+ channel [KN in P-wave). A single channel approach 
with enough time separation will be able to determine the ground state reliably. 
But if the pentaquark state is close to the KN scattering state, as is considered in 
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several calculations as a possbility in the 1/2 - channel, variational approach with 
multiple interpolation fields is a more effective approach. As a rule of thumb, the 
highest state will always be contaminated by still higher stat es and , thus, cannot be 
trusted. It is with the multi-operator variational calculations ^ ^ ^ ^ that one realizes 
that the earlier signal El of a low-lying l/2~ state in addition to the KN state at 
the threshold is spurious and is due to the fact that there were only two operators 
in the variational calculation. 

It is also widely speculated th at th e diqaurk-diqaurk-antiquark interpolation 
field proposed by Jaffe and Wilczek^O 

XJW = e^ c e de fe abc (u Ta C l5 d») 3, (u Td C l5 d & ) ltll5 Cs T ° (6) 

will have a better overlap with the pentaquark state to reflect the structure of the 
pentaquark. We note that this operator is very similar to the operator xt 

X A = e9 /c £d e/ £ a fc(/T75(it)(li ™^ 75(i e )>75C -T 9i (?) 

The only difference is that there is a covariant derivative between the two diquarks 
m X.JWi while X3 has a 7 M between u T and d in one of the u — d diquark pair. 

It turns out the pentaquark correlation functions built from these two inter- 
polation fields are linearly related. To see this, we consider the part of the cor- 
relator where the £> M and 7 M operate on the d quark propagator to the right 
which can be written in terms of the eigenstates of the Dirac operator, i.e. 
S(x,y) = J2 a ^ai^ip^iy) / (i\ a + Til). Now, the 2? M and 7 M operating on each eigen- 
state give 

The upper component oiT>jip a and "fjip a , i- e. T>jU a and —i(jjL a are related through 
the Dirac eigenvalue equation 

VjU a + it jkl a{D k U a = (Xj - W 4 )ajL a . (9) 

The lower components are similarly related. Thus, as far as the role of the interpo- 
lation fields is concerned and barring special exceptions, one would expect XJW to 
be equivalent to xt which is in term equivalent to xi an d Xi- There is no a priori 
reason why one operator will preferentially project to a particular state, be it the 
pentaquark or the KN scattering state. We should remark that the derivative op- 
erator xjw is motivated by the nonrelativistic picture and projects to the upper 
component of the Dirac eigenstates in the l/2 + channel; whereas, 75Xi;75X2> and 
X3 utilize the lower components. For a large quark mass, <7jL a is approximately 
equal to VjU a /m from Eq. @. In this case, the correlation function due to a local 
operator is 0(m 2 ) smaller than that of the derivative operator. This is the rea- 
son that the derivative operators are usually used to calculate the orbitally excited 
states in heavy quarkonia. To illustrate this point, we plot in Fig. 3 the spectral 
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Fig. 3. The spectral weights (in arbitrary units) of S\i and the Nrf ghost state as a function of 
rr4(GeV 2 ). 

weight of the 5*11 which is obtained from the lower component of the nucleon cor- 

H lil 

relation function with the 1 — 74 projection to the negative parity state We see 
that the spectral weight does go down like m 2 in the range of pion mass from 250 
MeV to 350 MeV. Thus, for heavy quarks, it is better to use the derivative operator 
as the interpolation field and obtain the negative-parity Sn state from the upper 
component of the Sn correlator. For small quark masses, the lower component of 
the nucleon correlator works for the Su state just as well. We believe the same is 
true with the pentaquark calculations. One shold be able to reach the same five- 
quark states with the derivative operator, the X3 operator, or equivalently the xi 
or X2 operators with comparable spectral weights when the quark masses are close 
to their physical values. 

3. Distinguishing Pentaquarks from KN Scattering States 

Since the 5-quark interpolation field will, in general, project to both the one- 
particle pentaquark states and the two-particle KN scattering states on the lattice, 
one needs to devise a way to distinguish them. To th is end , it was advocated to 
study the volume dependence of the spectral weig htsEDl. if it is a one particle 
state, the spectral weight for the correlator constructed with point source and zero 
momentum point sink has essentially no volume dependence. On the other hand, if 
it is a two-particle scattering state with relatively weak interaction, it is inversely 
proportional to the 3 volume from the normalization factor. The detailed study 
of such volume dependence with overlap fermions by comparing results from the 
16 3 x 28 and 12 3 x 28 lattices with pion mass m v as low as 180 MeV revealed 
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that the ground states of both the 1/2" and 1/2+ channels are scattering states^!. 
One needs to be careful here in view of the fact that it is more subtle to obtain 
the spectral weight than the mass of a state. The former is more sensitive to the 
fitting procedure and the length of time separation in the correlator. As it is shown 
recently that only when the time separation is large enough (~ 2.7 fm) and with 
accurate data is the volume scaling of the spectral weight for the 1 = 2 wir scattering 
state verified ^1 In the ifiV-pentaquark system that we are concerned with, the 
ground state mass is quite high compared to that of one hadron. Furthermore, the 
excited states of the system, which are the discrete KN scattering states with N 
and K each with discrete lattice momentum in units of 2tt / L, are stacked up more 
compact than the radial excitation in the one hadron case when the lattice size 
is reasonably large (e. g. more than ~ 2.4 fm). Both of these factors require a 
large time sepa rati on or a more sophisticated fitting routine, such as fitting with 
Bayesian priors El r with a variational approach, in order to obtain the spectral 
weights reliably El. The study of I = 2 hit in Ref. [14] serves as a caveat for the 
authors' earlier work on pentaquark where they found th e vo lume dependence of the 
spectral weight behaves like that of a one-particle state 

EH 

Another calculation El 

which observes little volume dependence of the spectral weight for the first excited 
state in the 1/2" channel, employed variational method with two operators and 
used two exponentials to fit the spectral weights. As is widely known in fitting 
procedures, one cannot trust the results of the highest fitted state which is the 
second state in this case. As we pointed out in Sec. [21 the lesson learned from the 
variational calculations with multi-operators^l vs two operators^ should be taken 
to heart. As such, the results of this study ^ should be taken with a grain of salt. 

Another way of distinguishing a scattering state from a pentaquark state is the 
clever idea of adopting the 'hybrid boundary condition' ^ with the anti-periodic 
spatial boundary condition for the u and d quarks and periodic condition for the 
strange quark. This way, the energies of the KN scattering states will be raised 
compared to those from the usual periodic boundary condition; while those of the 
uudds pentaquarks will stay the same^. Using this technique, it is found that the 
state near the KN threshold in the 1/2" channel is a scattering state. 

The third way to tell KN scattering state and pent aq uark state apart is to 
examine the volume dependence of the spectrum SSGEI Bl The one-particle state 
is not expected to be sensitive to the lattice volume when it is large enough for the 
specific quark mass; whereas, a weakly interacting two-particle state with relative 
momentum p will have a volume dependence since p is in units of 2ir/L. This 
approach applies well to our study since the KN interaction is weak and, as a 
result, it does not distort the discrete KN spectrum much from the non- interacting 
one. But it requires high statistics in order to discern the volume dependence which 
is milder compared to the 1/V scaling of the spectral weight. 
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4. Chiral Symmetry and Ghost States 

It is learned from studying the quenched chir al lo gs in m,, TOjv^l and the meson 
cloud effect in Roper resonance and Sn (1535) that chiral dynamics begins to 
play an important role in baryons when the pion mass is lower than ~ 300 MeV 
which characterizes the chiral regime. It is logically possible that the pentaquark 
state exists only for very light quark masses. After all, the first prediction of the 
anti-decuplet pentaquark is based on the chiral Skyrme model 0^1 Since most of the 
lattice calculations are carried out outside this chiral regime, the conclusion drawn 
from these calculations may not be relevant in this case. To have a definite answer 
on the existence of the pentaquark will require calculations in the chiral regime with 
dynamical fermions. When in the chiral regime, one needs to be concerned with the 
presence of ghost states for the quenched approximation and the partially quenched 
case when the sea quark and valence quark masses do not match. The relevant 
ghost state in our case is KNr} 1 with the would be rj loop not fully developed to 
that of a physical rj . Instead, it has a double 7r pole which has the pion mass and 
breaks unitarity by giving a negative contribution to the 5-quark correlator. The 
ghost state is a practical problem in the l/2 + channel in the chiral regime where 
the N, K and rj can be in relative S-wave which becomes lower than the P-wave 
KN state at certain mass and gives rise to a negative correlatorEl In this case, the 
ghost state and the physical states have to be fitted together before drawing any 
conclusion ED After the KNrj' ghost state is removed, it is found that the lowest 
physical state is the KN P-wave scattering state in the l/2 + channel El 

There is only one lattice calculation which claims to have detected a pentaquark 
in the l/2 + channel^. When extrapolated to the chiral limit, it approaches the KN 
threshold. This contradicts the findings in all the other calculations. We speculate 
that this could be due to the contamination of the KNr]' ghost state. Even though 
the quark masses are somewhat higher than the chiral regime, the fact that the 
spatial volume is small in this study (L = 1.8 fm) makes it more susceptible to the 
negative contribution from the the three-particle KNrj' ghost state which scales 
like 1/V 2 . Since the would be rj has the same mass as the pion, the ghost state is 
expected to reach the KN threshold at the chiral limit. Another noticeable feature 
of the calculation is that the overlap fermion results are found to be higher than 
those of the Wilson fermion for heavier than ~ 550 MeV in the positive parity 
channel and ~ 700 MeV in the negative parity channel^! In the 1 /2 + case, the mass 
difference can be as large as ~ 300 MeV. This needs an explanation. In contrast, 
ground state masses in the l/2~ channel from 6 different calculations with both the 
overlap and Wilson fermions are plotted in Fig. 28 in Ref. [8] and they pretty much 
lay on top of each other in the pion mass range from 420 MeV to 890 MeV. 

5. Conclusions 

To summarize, we have learned several lessons from studying the 5-quark sys- 
tem which includes two-particle scattering states in the spectrum. First of all, the 



1, 2008 5:11 WSPC/INSTRUCTION FILE Beijingpenta 



8 Keh-Fei Liu 

appearance of the interpolation field does not necessarily betray the structure of 
the hadron. To learn about the structure of a hadron, it is better to study the 
three- and four-point correlation functions. Secondly, it is essential to discern the 
one- or two-particle nature of the observed state before making a claim of having 
detected a pentaquark. In order to resolve states close by, it is necessary to perform 
variational calculation with more interpolation fields than the number of states of 
interest. Lastly, it is necessary to go down to small enough quark masses where the 
chiral dynamics dominates and get rid of the quenched or partially quenched ghost 
state in the l/2 + channel before one can be sure if the pentaquark exists. 

After close examination of the existing lattice calculations, we conclude that 
there is no convincing evidence that any pentaquark state has been observed so far. 

Of course, 'absence of evidence is not evidence of absence' To finally settle 
the issue, one needs to carry out realistic dynamical fermion calculations at physical 
quark masses, do a variational calculation to study the spectrum and the volume 
dependence of the spectral weight, and remove the potential ghost states in the 
partially quenched case. 
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